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In the physics of flavor mixing, the flavor states are given by superpositions of mass eigenstates. 
By using the occupation number to define a multiqubit space, the flavor states can be interpreted 
QQ , as multipartite mode-entangled states. By exploiting a suitable global measure of entanglement, 

■ based on the entropies related to all possible bipartitions of the system, we analyze the correlation 
' properties of such states in the instances of three- and four-flavor mixing. Depending on the mixing 
, parameters, and, in particular, on the values taken by the free phases, responsible for the CP- 

, ^ . violation, entanglement concentrates in certain bipartitions. We quantify in detail the amount 

' and the distribution of entanglement in the physically relevant cases of flavor mixing in quark and 

^ neutrino systems. By using the wave packet description for localized particles, we use the global 

. measure of entanglement, suitably adapted for the instance of multipartite mixed states, to analyze 

QQ ' the decoherence, induced by the free evolution dynamics, on the quantum correlations of stationary 

I neutrino beams. We define a decoherence length as the distance associated with the vanishing 

, of the coherent interference effects among massive neutrino states. We investigate the role of the 

I I ' CP- violating phase in the decoherence process. 

Ph| PACS numbers: 03.65.Ud; 12.15.Ff; 03.67.Mn; 14.60.Pq 
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G^. I. INTRODUCTION 

' Quantum entanglement as a physical resource plays a central role in quantum information and communication 

^ ] science As such, it has been mainly investigated in systems of condensed matter, atomic physics, and quantum 

^5 ■ optics. In fact, such systems offer the most promising possibilities of practical realizations and implementations of 

■ quantum information tasks. In the domain of particle physics, entanglement has been discussed mainly in relation 
i to two-body decay, annihilation, and creation processes, see for instance Refs. 0, |3, 0, d, 0, 0, B- In particular, 

attention has been focused on the entangled KqKq and BqBq states, produced in e+e" collisions [3, [l3|- Recently, 
the entanglement of neutrino pairs, produced in the tau lepton decay process t ^ Vr + -\- e~ -\- T>^ + Ve, has been 
analyzed in connection with the violation of Bell inequalities Q . 
' A fundamental phenomenon of elementary particles is that of particle mixing which appears in several instances: 
quarks, neutrinos, and the neutral if-meson system (lll.[l^. Particle mixing, consisting in a mismatch between flavor 
and mass, is at the basis of important effects as neutrino oscillations and CP violation [1J| . Flavor mixing for the case 
of three generations is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the quark instance [l^. 
?H ' and by the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) in the lepton instance Jij, _17]. The matrix elements represent 
the transition probabilities from one lepton (quark) to another. For example, the neutrino mixing is described by the 
following relation: 

Uel Ue2 Ue3 

U^l U^2 U^3 I I W2) I , (1) 

UtI Ur2 Ur3 

where the states jj/c) with a = e,fj,,T denote the neutrino flavor states, the states with i = 1,2,3 denote the 
neutrino mass eigenstates (with masses m^), and Ua,i denote the probability amplitudes of transition of the MNSP 
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matrix jji^'^'^SP)^ Analogously, for the quark mixing the CKM matrix connects the weak interaction eigenstates 
{\d') , \s') , \b'))'^ with the strong interaction eigenstates of the quarks |s), similarly to Eq. it results 

(|c?'), |s'), |&'))"^ — U^^^^^''{\d), \s), \b))'^ . From Eq. ([1]), we see that each flavor state is given by a superposition of 
mass eigenstates, i.e. \va) = + Ua2W2) + f^asli^a)- Let us recall that both and are orthonor- 

mal, i.e. (i^ali^p) — ^a,0 and = Sij. Therefore, one can interpret the label i as denoting a quantum mode, 

and can legitimately establish the following correspondence with three-qubit states: = |l)i|0)2|0)3 = |100), 

\i'2) = |0)i|l)2|0)3 = |010), ji^a) = |0)i|0)2|l)3 = |001), where \)i denotes states in the Hilbert space for neutrinos with 
mass rrii. Then, the occupation number allows to interpret the flavor states as constituted by entangled superpositions 
of the mass eigenstates. Quantum entanglement emerges as a direct consequence of the superposition principle. Let us 
remark that the Fock space associated with the neutrino mass eigenstates is physically well defined. Indeed, at least in 
principle, the mass eigenstates can be produced or detected in experiments performing extremely precise kinematical 
measurements. For instance, as pointed out by Kayser in Ref. [isj . in the process of pion decay 7r+ — > /i+ + I'fj,, 
highly precise measurements of the momenta of the pion and muon will determine the mass squared of the neutrino 
ml with an error Am^ less than the mass difference \mf — m|| {i ^ j = 1, 2, 3). Thus, the "physical neutrino" jz^i) 
involved in each event of the process is fully determined [l8l | . This kind of experiment will lead to the destruction of 
the oscillation phenomenon. 

Therefore, entanglement is established among field modes, although the quantum mechanical state is a single-particle 
one. This is in complete analogy to the mode entanglement defined for single-photon states of the radiation field 
or the mode entanglement introduced for systems of identical particles [l9i] : In all these instances, entanglement is 
established not between particles, but rather between field modes. In the particle physics instance, the multipartite 
flavor states can be seen as a generalized class of W states. The latter are multipartite entangled states that occur 
in a variety of diverse physical systems and can be engineered even with pure quantum optical elements (20| . From 
a theoretical standpoint, the concept of mode entanglement in single-particle states has been widely discussed in the 
literature and is by now well established [l^ [2l|, [S^, [23l. and linear optical scheme has been proposed to demonstrate 
multipartite entanglement of single-photon W states [24| . Experimental realizations include the teleportation of a 
single-particle entangled qubit [2&|, the quantum state reconstruction of single-photon entangled Fock states [26j . 
and the homodyne tomography characterization of dual-mode optical qubits using a single photon delocalized over 
two optical modes [2^ . Among the experimental proposals, we should mention a scheme for quantum cryptography 
using single-particle entanglement [28^. Moreover, remarkabl y, t he nonlocality of single-photon states has heen exper- 
imentally demonstrated by double homodyne measurements [29| , thus verifying a long-standing theoretical prediction 
[sol . [Slj . Very recently, the existing schemes to probe nonlocality in single-particle states have been generalized to 
include massive particles of arbitrary type [3^ , thus paving the way to the study of single-particle entanglement in a 
variety of diverse systems including atoms, molecules, nuclei, and elementary particles. 

Concerning the neutrino system, the main difference between the single-photon states and the single-particle neu- 
trino states is related to the spatial separability of modes. For instance, the polarization modes of polarization- 
entangled single-photon states can be easily spatially separated by means of a polarizing beam splitter. On the 
contrary, at present, it is not available a beam splitter analog for neutrinos. However, it is worth noting to recall that 
spatial separability and nonlocality are not necessary requirements for entanglement ^sj. Nevertheless, the spatial 
separation between massive neutrino states emerges in the dynamics of the free evolution in the wave packet approach 
[ssl . [U [ssl . [36| . In quantum theory localized particles are described by wave packets; moreover, during the free 
propagation, the different mass eigenstates {vi) in the packet travel at different speeds. Thus, the evolution leads to a 
spatial separation along the propagation direction (time delay) of the mass eigenstates and the difference between 
their arrival times at a given detector is observable [13, [ll] . The "decoherence" induced by the free evolution leads 
to a degradation and even to a complete destruction of the oscillation phenomenon [s^, [Sa, [s^ . It is interesting to 
investigate the influence of the decoherence on the quantum correlation of the multipartite entangled mass eigenstates. 

The issue of mode entanglement in single-particle states of elementary particle physics has been recently addressed 
by the study of the dynamical behavior of bipartite and multipartite flavor entanglement in neutrino oscillation [soj . 
In the present paper we characterize the correlation properties of the multipartite single-particle states that emerge 
in the context of lepton or quark mixing. These states turn out to be generalized W^-like entangled states. By 
resorting to a suitable measure of global entanglement, we analyze in detail their properties for different occurrences 
of flavor mixing and particle oscillations both in the quark and in the leptonic sectors. Furthermore, by using the wave 
packet description for the free-propagating neutrino states, we analyze the dynamical behavior of the multipartite 
entanglement in the phenomenon of neutrino oscillations. The paper is organized as follows: In Section [TT] we discuss 
the main aspects of different measures of multipartite entanglement. Following the approach of Ref. [40[ , we introduce 
a characterization of multipartite quantum correlations based on suitable entanglement measures for all the possible 
bipartitions of the A^-partite system. In Section IIIII we recall the formalism of flavor mixing in order to deflne 
generalized classes of three-partite VF^'^-' and four-partite PF'^'*^ states. In Section ITVl we study in detail the behavior 
of entanglement for the and W'^'^^ states as a function of the free phases in the case of maximal mixing. In 
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Section |Vl we apply the general formalism developed in the previous Sections to the quantification of multipartite 
entanglement in the most relevant cases of quark and neutrino flavor mixing. Finally, in Section IVTl bv using the wave 
packet treatment, we analyze the effect of propagation- induced decoherence on multipartite entanglement. 



II. MULTIPARTITE ENTANGLEMENT 



In this Section we will briefly discuss the problem of quantifying multipartite entanglement in relation to global 
aspects and statistical properties, and introduce measures particularly suitable for our purposes. For recent, detailed 
reviews on the qualification, quantification, and applications of entanglement, see Refs. [4lll42l. l43l| . Concerning bipar- 
tite pure states, entanglement is very well characterized by proper and efficient measures. In fact, for a bipartite pure 
state pi2 the von Neumann entropy E^pf = —Tri[pi log2 pi], for the reduced density matrix pi = Tr2[pi2], completely 



determines the amount of entanglement 4j| . For a given bipartition, Ey]\r has its maximum log2 d, where d denotes 
the minimum between the dimensions of the two parties. For bipartite mixed states, several entanglement measures 
have been proposed [45l. l46l. [4711 . Although providing interesting operative definitions, the entanglement of formation 
and of distillation f45| are very hard to compute. A celebrated result is the Wootters formula for the entanglement of 
formation for two-qubit mixed states [i^. An alternative measure, closely related to the entanglement of formation, 
is the concurrence (the entanglement of formation is a monotonically increasing function of the concurrence) [49| . The 
same difficulties of computation are encountered with the relative entropy of entanglement [46j| . At present a com- 
putable entanglement monotone is the logarithmic negativity E_\f, based on the requirement of positivity of the density 
operator under partial transposition Ej\f = log2 || pi2 ||i, where || • |ji denotes the trace norm, i.e. II O \\i=Tr[V0W] 
for any Hermitian operator O [4^ . The so-called bona fide density matrix pi2 is obtained through the partial trans- 
position with respect to one mode, say mode 2, of pi2, i.e. pi2 = Pi2^- Given an arbitrary orthonormal product basis 
Ki 1^2), the matrix elements of pi2 are determined by the relation (ii ,i2\pi2\ki , ^2) = (ii ,^2|pi2|fci ,32)- Obviously, 
for pure states such a measure provides the same results as the von Neumann entropy. 

The challenge of quantifying entanglement becomes much harder in multipartite systems. Important achievements 
have been reached in understanding the different ways in which multipartite systems can be entangled. The intrinsic 
nonlocal character of entanglement imposes its invariance under any local quantum operations; therefore, equivalence 
classes of entangled states can be defined through the group of reversible stochastic local quantum operations assisted 
by classical communication (SLOCCs) Such an approach allows to demonstrate that three and four qubits can 
be entangled, respectively, in two and nine inequivalent ways [5l|, [111- In particular, all three-qubit entangled states 
are related to two fundamental classes of states: the GHZ state \GHZ^'^^) and the W state IM^^'^-') [5l|, |53|. In the 
A^-partite instance, such states are defined as: 

IGFZW) = ^(|000...0) + 1111...1)), (2) 
v2 

|M^(^)) = ^(|100...0) + |010...0) + |001...0) + ...|000...1)). (3) 

V N 



The GHZ and W states are fully symmetric, i.e. invariant under the exchange of any two qubits, and greatly differ 
each other in their correlations properties. The GHZ state possesses maximal A^-partite entanglement, i.e. it violates 
Bell inequality maximally; on the other hand, it lacks bipartite entanglement. For instance, in the case A^ = 3, 
abandoning one mode, the resulting mixed two-mode state turns out to be separable. The W states possess less 
A^-partite entanglement, but maximal fT-partite entanglement {K < N) in the AT-reduced states. 

Several attempts have been done to introduce efficient entanglement measures for multipartite systems. The 
characterization of the quantum correlations through a measure embodying a collective property of the system, 
should be based on the introduction of quantities invariant under local transformations. A successful step in this 
direction has been put forward by Coffman, Kundu, and Wootters. Studying the distributed entanglement in systems 
of three qubits, they defined the so-called residual, genuine tripartite entanglement, or 3-tangle, a quantity constructed 
in terms of the squared concurrences associated with the global three qubit state and the reduced two-qubit states 
[49| . While successfully detecting the genuine tripartite entanglement in the state IGHZ^^")), the 3-tangle vanishes 
if computed for the state |VF^'^^), thus being not appropriate for this class of states. Several generalizations of the 
3-tangle have been proposed [s^ . The Schmidt measure, defined as the minimum of log2 r with r being the minimum 
of the number of terms in an expansion of the state in product basis, has been proposed by Eisert and Briegel [5^. The 
measure vanishes if and only if the state is fully product, thus it does not discriminate between genuine multipartite 
and bipartite entanglement. However, the Schmidt measure is able to distinguish the GHZ and the W states; for 
instance, it yields the value 1 for \GHZ'--^^) and the values log2 A^ for the |Vl^^^-') state (considering A^-partitions of 
the system) . Multipartite entanglement can be characterized also by the distance of the entangled state to the nearest 
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separable state; this is the geometric measure [56|. Simpler proposals are given in terms of functions of bipartite 
entanglement measures [13, [13, [H, [H, [gO] ■ An example of this kind of proposals is the global entanglement of Meyer 
and Wallach, that is defined as the sum of concurrences between one qubit and all others [13], and can be expressed as 
the average subsystem linear entropy [SS^ . A generalization of the global entanglement has been introduced by Rigolin 
et at, using the set of the mean linear entropies of all possible bipartitions of the whole system Hc?]. Recently, another 
approach has been proposed, based on the distribution of the purity of a subsystem over all possible bipartitions of 
the total system [60|- 



A. Average von Neumann entropy 

We intend to analyze the entanglement properties of a generalized class of W states (finite-dimensional pure states) . 
To this aim, we adopt an approach similar to that of Refs. [13, 113, [H, [El, [fi^l , thus characterizing the entanglement 
through measures defined on the possible bipartitions of the system. As we are dealing with pure states, we define as 
proper measure of multipartite entanglement a functional of the von Neumann entropy averaged on a given bipartition 
of the system. Let p = \ip){'ip\ be the density operator corresponding to a pure state {tp), describing the system S 
partitioned into N parties. Let us consider the bipartition of the A^-partite system S — {Si, S2, ■ ■ ■ , S^} in two sub- 
systems Sa„ = {5*1^, S'i2, . . . , S'i„}, with l<zi<«2<---<«n<A^(l<"-<A^), and S'_b„_„ = {Sj^,Sj^, . . . , S'-,„_„}, 
with 1 < ji < ^2 < • ■ • < jN-n < N, and ig 7^ jp. We denote by 

PA„ = P^,,^^,...,^„ = Trs„_„ [p] = T'^'ji , • • • JJV-„ [p] (4) 

the density matrix reduced with respect to the subsystem S'bjv-„- The von Neumann entropy associated with such a 
bipartition will be given by 

^(A„;S„_„) ^ j^^^ j^g^ ^^^j 

At last, we define the average von Neumann entropy 



where the sum is intended over all the possible bipartitions of the system in two subsystems each with n and N — n 
elements {1 < n < N). 

For instance, in the simple cases of a three qubit states, as the states pw(3) = | W^'^^) (W^^'^-' | and Pghz(^) — 
\GHZ^'^'>)(GHZ'^^^\, only unbalanced bipartitions of two subsystems can be considered. Straightforward calculations 
yield 

Eif EE £;^^^'^^)(p^^<3,) = {Ei%'\pwis,)) = log23- ^ ~ 0.918296, (7) 

Eiif^''^^\pGHZ(3)) = {eI%^\pghz(3))) = 1- (8) 

On the other hand, for a four-qubit state we have both unbalanced, i.e. {Sas, Sbi), and balanced bipartitions, i.e. 
{Sa2,Sb2)- For the state pyyw = \W(^^){W'^'^'>\, we get 

^ 4^^^^^) (pv^(4, ) = {Elf^ {p^,r, )) = 2 - J log2 3 ^ 0.811278 , (9) 

Eg^ ^ Eii;--''^\p^,.,) = {Eif\p^,.,)) = 1. (10) 

Of course, all the measures evaluated on the state Pghzw give the maximal, normalized value 1. It is worth noting 
that in order to characterize the multipartite entanglement in a iV-partite system, the number of bipartite measures 
grows with N. 



B. Average logarithmic negativity 

As well known, the entropic measures cannot be used to quantify the entanglement of mixed states. In order to 
measure the multipartite entanglement of mixed states, and following the same procedure as in the previous subsection, 
we introduce a generalized version of the logarithmic negativity [47[ . Let p be a multipartite mixed state associated 
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with a system S, partitioned into N parties. Again we consider the bipartition of the A^-partite system S into two 
subsystems Sa„ and S'b„_,^. We denote by 

=p^^^"-" = pPTn,n,...,j.^^ (11) 

the bona fide density matrix, obtained by the partial transposition of p with respect to the parties belonging to the 
subsystem S'b„_„. The logarithmic negativity associated with the fixed bipartition will be given by 

£;CA„;B„_„) ^ log2||p^Jog2P^„ 111 . (12) 

Finally, we define the average logarithmic negativity 



E<"^"-'^ (13) 



where the sum is intended over all the possible bipartitions of the system. 



III. GENERALIZED W STATES IN FLAVOR MIXING 



In this Section, we consider generalized W states of the form: 

N N N 

|WW(ai,a2,...,aAr)) - J] afc |<5i,fe, <52,fc, • ■ • , 5jv,fc) = ^ |j/f , ^|a,|2^1, (14) 

k=l k=l k=l 

where Sm,n denotes the Kronecker delta. In particular, we will consider the cases corresponding to A^ = 3, 4. Moreover, 
we will adopt a parametrization for {a^} commonly used in the domain of elementary particle physics, and associated 
with the phenomena of A^-flavor mixing, i.e. quark and neutrino mixing [ll| . 

The orthonormal set of flavor states |V';^^) are defined through the application of the N x N mixing matrix [/(^/) 
to the basis vectors li'^.^^), i.e. — J2k=i ^i'k^\^k^'^) {l,k — I, . . . , N). An N x N unitary matrix contains, in 

general, A^^ independent parameters. Each of the 2N fields (two for each lepton generation) can absorb one phase. 
Moreover, there is an unobservable overall phase, so we are left with (A^ — 1)^ independent real parameters. Among 
these, ^^^2 rotation angles, or mixing angles, and the remaining — — ^"^^^ are phases, which are responsible 

for CP violation. Applying this formalism, we determine A^ orthonormal flavor states {tpi^^), that belong to the class 
of generalized W states defined by Eq. p4)) . 



A. Generalized W'^' states from three-flavor mixing matrix 

In the case of mixing among three generations (either leptons or quarks), the standard parametrization of a 3 x 3 
unitary mixing matrix is given by [ll|: 



U{0,S)\iiJ (15) 




C12C13 S12C13 

U{9,S) = I -S12C23 - Ci2S23Sl3e*'' C12C23 - Si2S23Sl3e*'^ S23C13 I , (16) 
S12S23 ~ Ci2C23Sl3e*'' -C12S23 - Sl2C23Sl3e"' 

where \iif) = (|i^e), Wn), Wt})'^ are the states with definite flavor and |£^) = (|i^i), |j^2), Ws))'^ those with definite 
masses. In Eqs. (fT5|) and (fT6|) . the following shorthand notation has been adopted: (0,6) = (6*12, 6'i3, 6*23; (5), = 
cos 9ij and Sy = sin^y . In this case, we have three mixing angles 9i2, ^13, ^23, and a free phase i5. It can be shown 
that the values of these parameters for which the three flavor mixing is maximal are [l5| : 

nrnax _ . nmax _ . nmax _ „_.„„„„ , / ^. cmax ^^^\ 

9i2 - ^23 - ^) ^13 -arccosW-, d --. (17) 
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In correspondence of these values, the matrix elements in Eq. (|16p have all the same modulus 

For = 3, we define the generalized class of three-qubit W states as those generated by means of the following 
matrix, which is obtained by the above mixing matrix upon multiplication of the third column by e^^: 

iWj'^Ho-.S)) = C/(3/)(^,^)|iy(3)) (18) 

/I \ 

t/(3/)(^,5) = u{e,s) 1 , (19) 




where |W:(3)(^;^)) = (| W^i^^^, 5)), | W^i^^ (^, 5)), I W^i^^ (0, ^))) ^ and = (kf ^), '))^. The entangle- 

ment properties of the states associated with matrices (IT51) and are identical. When all the mixing parameters 
are chosen to be maximal as in Eq. (jl7p . the matrix U^^-^^ becomes: 

U^ll I . (20) 

with y = exp (2i7r/3). In the case of maximal mixing, all the states possess the same entanglement of 

£^i^'''^'^(lii:^^n6'""'^'5'""^))) = (^;f^,^^(|]£^^H^™''"^;5'"°^)))) = E^^ , (21) 

where £^21'' is defined in Eq. ([7]). In the next Section we will analyze the entanglement properties of the \wl^\6,5)) 
states, and their behavior as a function of the mixing parameters. 



B. Generalized W'*^ states from four-flavor mixing matrix 

Let us now consider the four-flavor mixing [N = 4). In particle physics, such a case could be realized, for instance, 
by a situation in which there are three active neutrino types and an extra one, non interacting, the so-called "sterile" 
neutrino. Obviously, such states correspond to physically realizable situations in optical and condensed matter 
systems. The corresponding four-flavor mixing matrix U'^'^^\9^ 5) will be built on 9 independent parameters, 6 mixing 
angles and 3 phases, i.e. {9; 6) ~ (^12, ^13, ^14, ^237 ^24j ^34; <^147 <^237 <^34)- Explicitly, the mixing matrix for four flavors 
can be written as the following product of elementary matrices: 



u^^f\e-~5) = £734(^34, ^534) £724(^24)^^23(^237 <^23)C7i4(0i4, 5u)Ui'i{0i'i)Ui2{9i2)'Us{5i4) , 



(22) 



where 



Us{5ia) 



/ 1 

10 

1 

Vo e^'^i" 



Ur 



cos 612 sin 6*12 
— sin 6*12 cos 6*12 



1 
1 



Ui 



cos6'i3 sin 6*13 

10 

— sin 013 cos(?i3 

1 



(23) 



Ui. 



V 



cos 6*14 e"**" sin 6*14 ^ 
10 

1 



23 



sin cos 014 J 



/I 








COS 023 
j'^^s sin 6»23 









6"^''=^ sin 023 



cos 023 







1/ 



(24) 



U24 



/I 

COS 024 sin 024 

10 

\ — sin 024 cos 024 



34 



/I 

1 



\0 






cos 034 e 
gi<534 sin 6)34 







-'^^'^ sin 034 

cos 034 



(25) 



Analogously to the definition (jlSp given in subsection llll Al the class of generalized four-qubit W states can be defined 



as 



(26) 
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The matrix (|22|) is maximal, i.e. all elements have the same modulus 1/2, for the following set of values: 



9 max nniax 
12 — "34 



(>14 = <P, 



vr 
4' 



9 max nraax 
14 — '^23 — 



9 max 
13 



arccos \j - ; 



9 max 
24 



arcsm ^/ -; 



cmax _ 
(^23 = TT - 

(4/), 



(^34 = 



For the choices ([27| and ([SS]) . Umax{4>) takes the simple form 



/II 1 1 



1-11 1 
\ 1 -1 -e*"* e"^ 



All the states \W^*\9 



(27) 
(28) 

(29) 



max . j:max 



)) exhibit the same amount of entanglement of the standard four-qubit jW^'-'*^) state: 



E, 



E, 



(A2\B2) 

'N 



Qyy{i)(^Qmax.Jmax^^^ = (i;^^^^^ ( | W^^"*) (fii""^ ; )) )) =E^fK 

(|W^('')(^""^;(5'""^))) = (£;^^^^(|VK(4)(6i""'"=;5""^)))) = i;^2^ 



(30) 
(31) 



for any bipartition {A2 ,-62) and (A3 ,i3i). E^'^^ and £'22'' are given in Eqs. ^ and pU|) . respectively. 



(4) 



IV. THE CORRELATION PROPERTIES OF jW^^^) STATES 



In this Section we analyze the correlation properties of the class of W-like states defined by Eqs. (fT8| and (|26p . Such 
properties are completely determined by the free parameters of the mixing matrix formalism, i.e. the rotation angles 
9ij and the phases 6ij. Let us recall that for iV = 3 we have three angles and one phase, while for iV = 4 we have 
six angles and three phases. In our formalism, the state associated with the first row of the matrices U^^^\9,S) and 

U^'^^\9,5), i.e. the states \We^\9\5)) and \We'^\9]5)), respectively, reduce to standard 3-qubit and 4-qubit W states 
by fixing the rotation angles to their maximal values 9^°"^, according to Eqs. ^T7\ and (|?7|) . Therefore, in the instance 

of A^-fiavor W states with maximal mixing angles, i.e. \W_^-'^\9"^'^^; S)), there exists a subspace (of dimension N — 1), 
that is is orthogonal to the state and is spanned by the vectors ^(fi*™"^; ^)), . . . |M^i^^(6i™"^; ^))|. For 

simplicity, in the following we will restrict ourselves to the study of the entanglement properties of such a subclass of 
generalized states, which are parameterized by the phases of the mixing matrix. 



A. Case of \W}^\e""''' ■ S)) states 



First, we discuss the entanglement properties of 3-partite W states U''^^\9"^°'^ , 5); in particular, we study the 
dependence of entanglement on the phase 6, with the rotation angles 9ij at their maximal values 9Yj°'^ , given by 

Eq. p7|) . In this way, we obtain a set of three orthogonal generalized W states \Wa {5)) = \Wa (6*™°^; (5)) (a = 
e,iJ,,T), of which the first one is the usual |IF^^^) state. Correspondingly, the matrix U'^'^^^ is specialized to 



1 

7! 



Let us compute the quantities E^^'^^'^ and (iJ^ 



E, 



/I I 1\ 

-i(\/3 + e**) i(V3-e'*) e'^ . 
V i(\/3-e*^) -i(V3 + e'*) e^V 

as defined by Eqs. (O and ^ in Section Hi Al We get: 



(2:1) 
N 



i(l:2;3) 

vNe 


p(l:3;2) 
~ ^vNe 


p(2:3;l) 
~ ^vNe 


p(l:2;3) 


n(l,3;2) 
vN (1 


_ p(2,3;l) 




cos6\ , 
2V3j^°^ 


1(2, 3;1) 
vN ^1 


_ r(1,3;2) 




cos(5\ , 
2^/3j'°^ 



= E. 



(1,2;3) 
vN T 



= log2 3 - - 



"1 


cos 6 




cos (5 


3 ^ 


2^3. 




2V3 


"1 


cos b 




cos (5 


3 + 


2\/3. 




2\/3 



l0g2 
l0g2 



COS (5 



2V3 
cos 5 

2V3 



(32) 



(33) 
(34) 
(35) 
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FIG. 1: (Color online) The von Neumann entropy E^^'^'^'' and the average von Neumann entropy (S^^^^) as functions of the 
CP-violating phase 5. E^''^'^^ is plotted for the following bipartitions k: (a) i = 1, j — 2, and fc = 3 (dotted line); (b) i — 1, 

is constant and takes the reference value 



3, and k = 2 (dashed line); (c) i = 2, j 



3, and k = 1 (dot-dashed line). E'^^'^^ 



E^f = 0.918296. The average entropy {E^j^ ') (full line) attains the maximal value E21 at 5 = ipyr, with p integer 



p(2:l) 



(3) 



where the superscript {i,j;k) explicitly represents the specific composition of the bipartitions A2 = {Si,Sj} and 
Bi = {Sk}, with i,j,k — 1,2,3 and i ^ j ^ k. Moreover, in order to simplify the notation, the definition e'^^'^'^ = 

E'^^N^Hlwi^^S))) has been introduced. The states \wi^\S)), with a — fj,,T, possess correlation properties dependent 
on 6. 

Let us, for instance, consider the state |W^^^^((5)); in Fig. [U the plots display the behavior of e'^^"^'' and {eI^j^J^} 
as a function of S in the range [0, 2tt]. While E)^^'^ ' (dotted line) takes the constant reference value i?2i (as for the 
state IM^*^"^^)), the quantities E^j^'^'^ (dashed line) and E^^j^'^^ (dot-dashed line) vary with 5, attaining the absolute 



maximum 1 at the points (5i = ± arccos y'^j ^ 2p7r and (52 = ± arccos ^ 2p-n (with p integer), respectively. 

Therefore, the state jW^i {^i))^ with i = 1,2, exhibits maximal entanglement in a given bipartition, equal to the 
entanglement shown by the GHZ state \GHZ''^^). Moreover, for each given range of values of 5, we see that either 

(dot-dashed line) exceeds the reference value E'f^ . This phenomenon of periodic 



e: 



(1,3;2) 



p(2,3;l) 



vNt. (dashed fine) or E\^j^ ^ 



entanglement concentration is reminiscent of spin squeezing in collective atomic variables. On the other hand, the 
average von Neumann entropy {E^^j^ p stays below the reference value i?2i j attaining it at the points ^ = f ±p7r. 
In conclusion, the free parameter 6 can be used to concentrate and squeeze the entanglement in a specific bipartition, 
allowing a sharply peaked distribution of entanglement, corresponding to a lowering of the average von Neumann 
entropy. 



B. Case of \W}^'> {9"''''' ; 5)) states 



Due to the increased number of degrees of freedom, the class of VK-like states for = 4, i.e. Eq. ((26|) . yields a more 
complex scenario for investigation. Proceeding as in Section [IV Al we fix the rotation angles at their maximal values 
given by Eq. (P7)) . and leave free the phases Sij. The matrix U^^^^{6; S) acquires the form: 



qmax 



/ 



v 



1 

1. I ^23 ^14^34 i , Z34 ]_ _ 

2~^2 ^3'^6~^2 2 

]. ^14 I -^23 2:232:34 _|_ 2:34 I 

2 3^6 2~'~2 2 



1 

-223 Zl4Z2_^ . ^23 ^3-^ 

2 § 6 

— ^14 I ^23 I ^23^34 
3 6 ~'~ 2 



Z34 

2 



_ ^14 I ^^23 
214-234 -223-^34 

3 . 3 

^34 



Zl4 

3 



^23 



1 

Zi4 

^34 •^14 
Zi4 



where 



(36) 



The explicit analytical expressions for the entanglement measures evaluated on the states 
|Wi'*''((5)) = \wi^\9"^'^^ ; S)) {a — e,^,T,s) are rather long and involved, and are reported in Appendix El Note 
that the state |VFi''''((5)) coincides with the usual jVK'^'^') state. As an example, let us analyze in detail the entangle- 
ment of the state \wlf\6)), that depends on the phases ^14 and S23 and is independent of the phase ^13. In Fig. [21 
the plots I-III display eI]^'^'^\ and e'^]^^'^'^\ respectively, as a function of and 623] the plot IV displays 
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(i i'k I) 

FIG. 2: (Color online) The von Neumann entropy E^^'^' for balanced bipartitions and the average von Neumann entropy 

{^vN fi) functions of the phases S14 and 523- Panel I shows S.),^^ .It exhibits an osciUating behavior along the direction 

parallel to the vector (5i4, 523) = (1, 1). Panels II and III show the entropies e'^^'^''^'^ and E^^m^'^'^\ respectively. They exhibit 
a nontrivial behavior yielding a periodic array structure of holes and dips. The combined behaviors of all the entropies result 
in the average von Neumann entropy, displayed in panel IV. All the four functions reach the maximum attainable value 1 of 
the entanglement at (5i4 + ^23 = ±P7r, with p odd integer. 



the behavior of the average entropy (-E.^jy^). The entanglement takes the maximum value 1 in correspondence of the 

values given in Eq. (p8|) . i.e. for (5i4 + 523 — with p odd integer. Moreover, while -E^jv^^f exhibits an oscillating 

behavior along the direction parallel to the vector ((5i4,(523) = (1, 1), the quantities E^P4"I^'^\ El]f'^'^\ and (eI^j^J^) 
show a periodic array structure of holes. 

fr 7 k I) 

Next, we consider the entropies corresponding to the unbalanced bipartitions E]^^'^' . The surface plots of these en- 
tropic measures, as functions of ^14 and ^23, are similar to those for the case of balanced bipartitions, shown in Fig. 
In order to better highlight their structure, in Fig. [3l we plot one-dimensional sections of the surfaces belonging to 
the plane 5u = (523- We see that, as in the three qubit instance, concentrations of entanglement (with a value in the 
range [e'^ , 1]) occurs for the bipartitions (1; 2, 3, 4) and (2; 1, 3, 4), corresponding to a lowering of the average entropy 
{E^^\). In the range [0, 27r], both e'^^'^''^^ (dotted line) and E^j^'^'^^ (dashed line) exceed in alternating order the 
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FIG. 3: (Color online) One-dimensional sections of the von Neumann entropies -B^'j^^'*' for unbalanced 1 : 3 bipartitions and 

their average as functions of the phase 5 (5 = 5i4 = 523). Similarly to the case of three qubits, in the unbalanced four- 

qubit instance a concentration of entanglement can be observed in the entropies B'^jv^^'^^' (dotted line) and E^^^'^''^^ (dashed 

line). The entropy E''^^'^'^^ (double-dot-dashed line) is constant and takes the reference value E'^ = 0.811278. The entropy 

-^uiv^/i^'"'' (dot-dashed line), and the average entropy {El]^j^) (full line) are always limited by this value, reaching it at points 
(5=1 +pn. 



reference value E^^^ , and attain their maximum value 1, respectively at the points Sa = ± arccos §(\/2 — 1) 



± 2p7r 



and dh 



± arccos 



|(V2 — 1) ± 2p7r. This behavior is again reminiscent of spin squeezing in atomic systems. 
Analogously to the three-qubit instance, the average entropy exhibits an oscillatory behavior, and stays below the 



.(4), 



reference £^31'', reaching it at ^ = ^ -I- p7r. 

As they depend non trivially on all the phases Sij, the states |iyr^''((5)) and |Wi''''((5)) possess an even richer structure 
of quantum correlations, compared to the case |M^/1^^(5)) However, in both instances, one observes similar effects as 
the ones that occur for the state \wlf \S)). 



V. QUANTIFYING ENTANGLEMENT IN QUARK AND NEUTRINO FLAVOR MIXING 



In this Section, we quantify the entanglement in situations of quarks or neutrino mixing, described by the three 
flavor states defined in Eg. (f T5 |l . We will set the parameters of the matrix (fT6|) at the values established by the current 
experiments [l^ [61I [6^. [63j . In the case of quarks, the mixing angles of the CKM matrix, are given by p^ : 



qCKM 

'12 



9g^^' = 0.2° ±0.1°, 



9g,^*^ = 2.4° ±0.1°. 



= 13.0° ±0.1°, 

Moreover, a measurement of the CP violation has yielded the value for the CP- violating phase [IS 



:CKM 



(37) 
(38) 

In Table H we list the values for the von Neumann entropies -E^jv a ^ 1 with a — d' , s', b' and i, j, k = d,s, b, and {E^^^P^ 
corresponding to the states (fT5|) . with the mixing angles and the CP-violating phase fixed to Eqs. ((37|) and (|38p. 
respectively, without taking into account the uncertainties. We see that, in the range of the experimentally measured 



= 1.05 ±0.24. 



Q 


p{d,s;b) 
'^vN a 


^vN a 


p(s,6;d) 
^vN a 


/p(2:l)\ 


d' 


0.0002 


0.2889 


0.2890 


0.1927 


s' 


0.0185 


0.2960 


0.2887 


0.2011 


b' 


0.0186 


0.0180 


0.0010 


0.0126 



TABLE I: von Neumann entropies El''j^'^^ and {E^^p^lj {a = d',s',b') for the three-flavor states associated with the quark 
mixing. 

values of the mixing angles, the entanglement stays low, very far from the maximum attainable value 1. Moreover, it 
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concentrates in the bipartitions (d, b; s) and (s, b; d) of the states \d') and |s'), while it is very small for the state \b'). 
In the case of neutrinos, the most recent estimates of the parameters of the MNSP matrix are expressed by the 
following relations [63l |: 



sm 



2 nMNSP 
'12 



0.314(1 



-0.18 
-0.15 



sm 



qMNSP 
^13 



= (0., 



-2.3 
-0.8 



X 10" 



sm 



qMNSP 
-'23 



0.45(1 ^'^•^^). (39) 



The CP-violating phase associated with lepton mixing is, at present, completely undetermined; therefore, J^^^sp 
may take an arbitrary value in the interval [0, 27r). In Table [III by using the relations (|39p (without taking into 
account the uncertainties) and for arbitrary S^^^^'^ ^ we list the entropies corresponding to the neutrino flavor states. 
The given intervals of possible values are obviously due to the freedom in the choice of the CP-violating phase. 
Comparing Tables [J and |TT1 it turns out that the neutrino mixing states are more entangled and their entanglement is 



Q 


p(l,2;3) 


p(i,3;2) 


p(2,3;l) 
^vNa 


/p{2:l)x 
\^vN a/ 


e 


0.0672 


0.8948 


0.9038 


0.5995 




0.9916 


0.9220 - 0.9813 


0.5679 - 0.7536 


0.8469 - 0.8891 


T 


0.9939 


0.8397 - 0.9352 


0.4784 - 0.6922 


0.8025 - 0.8419 



TABLE II; von Neumann entropies E^^'^ 
mixing. 



and {Ei%'l) (a 



t) for the three-flavor states associated with the neutrino 



more homogeneously distributed among the different bipartions, compared to the quark mixing states. In the case of 
neutrinos, the uncertainties are very large. Moreover, the value taken by the mixing angle O^^^^^ is crucial. In fact, 
only if such an angle is non-vanishing, then the entropies are dependent on the CP-violating phase. Therefore, it is 
interesting to investigate the behavior of entanglement when one takes into account the experimental uncertainties 
on the mixing angles. To this aim, we assume that OfJ^^^ takes random values normally distributed around the 

experimentally observed values. For instance, in Fig.[4l we plot eI^^'I^^ and {E^p^j^) as a function of the free parameter 
gMNSP ^ g -y^g -(-j^g^-f- |-j^g entanglement corresponding to bipartitions (1,2; 3) and (1,3; 2) keeps high, (panels I 
and II); on the other side, the bipartition (2,3; 1) exhibits lower amount of entanglement (panel III), leading to a 
lowering of the average amount of global entanglement (panel IV). Thus, we can conclude that, for the states ji'p), 
the parties 2 and 3 are more strongly correlated compared to the pairs 1,2, and 1,3. Similar conclusions hold for the 
states |i^e) and Ii^t)- 



VI. DECOHERENCE IN NEUTRINO OSCILLATIONS 



In the previous Sections, we have provided an analysis of the mixing effect in terms of the quantum correlations 
of multipartite mode-entangled states, by exploiting tools commonly used in the domain of quantum information 
theory. The characterization of the entanglement of generalized multipartite W states, through the measurement of 
the amount of the quantum information content of these states, constitutes a description of a fundamental effect of 
particle physics. The physical insight of such an analysis acquires even more relevance if it is transferred to a dynamical 
scenario, by studying the phenomenon of particle oscillations. Let us recall that both the phenomenon of particle 
oscillations and quantum entanglement are due to the superposition principle which gives place to coherent interference 
among the different mass eigenstates. In the particular instance of neutrinos, the standard theory of oscillations is 
developed in the plane- wave approximation [6J]. Adopting such an approximation, all the results obtained in the 
previous Sections hold for any time in the free evolution dynamics. However, a more realistic description of the 
phenomenon can be achieved by means of the wave packet approach [33, 34, 35, 36], for reviews see Refs. [37, 3^. The 
three massive neutrinos possess different masses; consequently, the corresponding wave packets propagate at different 
speeds, and acquire an increasing spatial separation with respect each other. Therefore, the free evolution leads to 
a natural lowering of the coherent interference effects, associated with the destruction of the oscillation phenomenon 
and with the vanishing of the multipartite quantum entanglement. In this Section, we intend to analyze the quantum 
correlations of multipartite entangled neutrino states by using the wave packet description for massive neutrinos. In 
particular, we want to study the "decoherence" effects, induced by the free evolution, on the multipartite entanglement 
among neutrino mass eigenstates. Let us notice that the forthcoming analysis, as well all the formalism developed in 
this work, can be applied to any system exhibiting the particle mixing. 

Following the procedure developed in Refs. [s^^fssl. [36|. by considering one only one space dimension, a neutrino 





FIG. 4: (Color online) The von Neumann entropies -E^'^J^'' for all possible bipartitions and the average von Neumann entropy 
{-^vN fi) functions of the CP- violating phase S. In panel I we plot the entropy E^j^ . It is constant and close to 1, the 
maximum attainable value of entanglement. In panel II we plot the entropy E^^^^I^K It is moderately (5-dependent and reaches 

its maximum at 5 — tv, still below 1 (within the experimental statistical errors). In panel III we plot the entropy E'^^^'^\ It 
corresponds to the bipartition with the least content of entanglement, is strongly dependent on 5, and reaches a minimum at 
(5 = TT. The resulting average entropy (E^^j^j^), displayed in panel IV, is weakly J-dependent and reaches a minimum at 5 = vr. 
The mixing angles 9^'^^^ are assumed to be Gaussian random variables, with a distribution centered at the mean values 

—MNSP 56^^^^^ 

9ij fixed to coincide with the experimental values (|39p . and a standard deviation aij chosen to coincide with — '-^ . The 

uncertainties &6f'-'^^^ are fixed at the maximum values between the left and right extrema given in Eq. (|39|l . The thick full 

T 1 • • ^ r\ A'J N ¥* "77-^-^ ^ S P 111 

lines represent the entropies with — 9ij , and null uncertainty. 



with definite flavor, propagating along the x direction, can be described by the state: 

W{x,t)) = Y,Uo..,^j{x,t)\v,), (40) 

3 

where the Ua,j denotes the corresponding element of the mixing matrix, is the mass eigenstate with mass m-,-, 
and ipj{x,t) is its wave function. Assuming for the momentum of the massive neutrino a Gaussian distribution 
ipj{p), the wave function is given by: 

^,{x,t) = ^ I dpV,(p)e^''-"»^^(^)*, ^P.ip) ^ ^^-A^e"^^"""^^', (41) 

where pj is the average momentum, ap is the momentum uncertainty, and Ej (p) = The density matrix 

associated with the pure state Eq. (^0]) writes: 

Pa{x,t) = \iya{xA)){iyaix,t)\. (42) 
If the inequality ap ^ E'j(pj)/mj holds, the energy Ej{p) can be approximated by Ej{p) ~ Ej + Vj{p — pj), with 
Ej = + m'j, and Vj = ^^q^^ \p-p- ~ ^ group velocity of the wave packet of the massive neutrino \vj). In 
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this case, the integration over p in Eq. (|4T|) is Gaussian and can be easily performed, yielding the following expression 
for Pa(x,t) 

Pa{x,t) ^ —^==2_^UajU^ke ^"^ W])Wk\, (43) 

where = (2o'p)^^- In the instance of extremely relativistic neutrinos, the following approximations are usually 
assumed 

9 9 2 

m~ rri; 

E.c^E + f^, p,- ~ (1-0— , f, 1 ^ (44) 

J ' ^ ^ ' 2E ^ 2EJ ^ ' 

where E is the neutrino energy in the limit of zero mass, and ^ is a dimcnsionless constant depending on the 
characteristic of the production process [H, [3^ . The density matrix (|43p provides a space-time description of neutrino 
dynamics. However, in realistic situations, it is convenient to consider the corresponding stationary process, which is 
associated with the time- independent density matrix pa{x) obtained by the time average of pa{x,t) [s^. By taking 
into account Eq. (|44p . and by computing a Gaussian integration over the time, the density matrix becomes [35| 



Pa{x) = ^C/ajC/*feexp 



2E 



ap 



W3){M, (45) 



with Am|j. = m'j — . The density matrix (j45p can be used to study, in the wave packet approach, the phenomenon 
of neutrino oscillations for stationary neutrino beams [s^, [H, [s^ . 

Here, we intend to analyze the coherence of the quantum superposition of the neutrino mass eigenstates, by look- 
ing at the spatial behavior of the multipartite entanglement of the state By establishing the identification 
Wi) = |<5i,i)i|(5i,2)2|<5i,3)3 = |(5i,i<5i,2<^i,3) (« = 1,2,3), we Can easily construct from Eq. (|45| the matrix with elements 
{lmn\pa{x)\ijk) , where i, j, k,l,m,n = 0,1. Let us notice that the density matrix pa{x) describes a mixed state, 
whose non-diagonal elements are suppressed by a Gaussian function of x. An appropriate quantifier of multipartite 
entanglement for the state pa {x) is based on the set of logarithmic negativities defined in subsection IIIBI We an- 
alytically compute the quantities E^^^''\ for i,i,k ~ 1,2,3 and i ^ j ^ and the average logarithmic negativity 
i^jVa^^ for the neutrino states with flavor a — e,fi,T. We assume for the mixing angles 9f^^^^ the experimental 



values l|39p. The squared mass differences are fixed at the experimental values reported in Ref. [63 



5m? -2 '9 



5m'^ = 7.92 X 10"^ eV"^ , 5m'^ = 2.6 x 10~^ eV"^ . (46) 

The parameters E and ap in Eq. (|^5)l are fixed at the values E = 10 GeV and ap = 1 GeV . Moreover, although 



depending on the particular production process 65| , the parameter ^ is put to zero for simplicity. In Fig. [5l we plot the 



logarithmic negativities for the electronic neutrino, i.e. -E^J-'g''^'' as function of the distance x. The bipartitions (1,3; 2) 
and (2,3;1), see panel I, exhibit a high entanglement content (> 0.93) that keeps almost constant for x < 10® m; 
finally, it goes to zero for a; ~ 3 x 10^ m. The bipartition (1, 2; 3) exhibits a low entanglement (< 0.24), that goes to 

zero for x w 9 x 10^ m. Furthermore, let us remark that the the logarithmic negativities and {E^^']}) for the 

electronic neutrino are independent of the CP-violating phase 5. 

In the muonic and tauonic instances, the independence from the CP-violating phase 5 holds no more. Therefore, 
first we choose to study the quantum correlations of these states for 5 = 0; then we consider separately the influence 
of a non-zero 5. In Fig. [51 we plot the logarithmic negativities for the muonic and tauonic neutrinos as functions of 
the distance x with 5 ~ Q. We see that the spatial behavior of multipartite entanglement for muonic and tauonic 

(\ 2'3) (\ 2'3l 

neutrinos are similar. The logarithmic negativities Ej^'^' and Ej^'^' are initially close to 1, and they go to zero for 

X ~ 10® TO. On the other side, -E'jJ^'^'^'' , E^'^'^'' , sjj^'^'^'' , and E^'^'^'' exhibit alternating regimes with slowly decreasing 
slope and with rapidly decreasing slope; moreover, all vanish for a; ~ 3 x 10® m. 

The average logarithmic negativity {Ej^'J) can be used to define a decoherence length Ldecoh as 

Ldecoh ■■ {E^^-^J) [Ldecoh) = 0. (47) 

From Figs. [51 [HI for assigned experimental parameters, we see that the common decoherence length for the neutrinos 
of flavor a — e, p,T can be estimated at a value of Ldecoh ~ 3 x lO^Km. 
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FIG. 5: (Color online) The logarithmic negativities -E'^'g''"' for all possible bipartitions and the average logarithmic negativity 

{E\f'J) as functions of the distance x. The quantities E\f\' (dashed line) and E\f\' (dot-dashed line), see panel I, show 
a high amount of entanglement content in the corresponding bipartitions, and seem to be superimposed. In panel II we plot 
a zoom of E]^'^' and E)^'^' to observe the differences in their behaviors: the two curves are initially separated, and then 
they superimpose each other. The bipartition (1,2; 3), associated with the quantity E^'^'^^ (dotted line), exhibits the lowest 
amount of entanglement. The full line corresponds to the average logarithmic negativity (sj/'P). The mixing angles O^J'^^^ 
and the squared mass differences Amfj are fixed at the experimental values (|39p and (|46|l . respectively. We assume the values 
E = lOGeV, Op = 1 GeV, and ^ = for the remaining parameters in Eq. (|45|l . All the plotted quantities are independent of 
the CP-violating phase S, that can be assumed arbitrary. The x axis is in logarithmic scale, and the dimensions are meters. 

Finally, we consider the influence of a non-vanishing phase S in determining the spatial behavior of multipartite 
entanglement of stationary neutrino beams. To this aim, in Fig. [7] we plot the logarithmic negativities for the muonic 
neutrino -E)/'^' and E)^'^' , with S fixed at the values 6 = 0, |, tt. The behavior of -Ej^'^' is not reported as it is 
independent of 6. We observe that the CP-violating phase S does not lead to a change of the decoherence length Ldecoh- 
However, we see that it may lead a lowering or an increasing of the amount of entanglement in a given bipartition, 
in agreement with the results obtained for the instance of static neutrinos. Similar results can be obtained for the 
tauonic instance. 



VII. CONCLUSIONS 

The study of entanglement between field modes can be fruitfully applied to a large variety of quantum mechanical 
systems, either in the usual case of many-particle multipartite entangled states or in the more intriguing instance 
of single-particle multipartite entangled states. In the present paper, stimulated by recent work on single-particle 
nonlocality and entanglement in quantum optical systems, we have extended the analysis of mode entanglement to 
systems of elementary particle physics. In particular, we have determined and studied the structural properties of the 
multipartite entangled states that occur in the physics of fiavor mixing, either in quark or in leptonic systems. These 
states are generalizations of the well known W states, endowed with nontrivial relative phases. These states include, 
as a special instance, the symmetric W state and the set of states orthogonal to it. We have implemented global 
and statistical approaches, based on the distribution of different bipartite entanglements, to quantify the generic 
aspects of multipartite entanglement in such states. We have studied in detail the correlation properties of three- and 
four-flavor W states. For properly chosen mixing parameters, we have shown that the phases, responsible for the CP- 
violation effects in particle physics, can be used to concentrate the entanglement in a particular bipartition, and we 
have identified some periodic patterns of entanglement concentration, dispersion, and revivals, that are reminiscent of 
spin-squeezing phenomena for the collective variables of many-body atomic systems. Moreover, we have analyzed the 
entanglement for the three-quark and three-neutrino mixing. In the particular instance of neutrino mixing, we have 
determined the effects of the free relative phases on the distribution of entanglement. By exploiting the wave packet 
treatment for neutrino mass eigenstates, we have considered in detail the influence of decoherence induced by the 
free evolution on the multipartite entanglement. A decoherence length can be defined as the distance associated with 
vanishing average global entanglement. Finally, we have studied the role of the CP-violating phase in the dynamics 
of free propagation. 
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FIG. 6: (Color online) The logarithmic negativities -E^'^'*' for all possible bipartitions and the average logarithmic negativity 
i-^N'a)' with a — (I, T, as functions of the distance x. In panel I we plot the negativities for the muonic neutrino. The 
bipartition (1,2; 3), associated with the quantity -E^J^'^'^' (dotted line), shows the highest initial amount of entanglement, 

that goes to zero for a lower of x with respect to the other bipartitions. Ej^'^' (dashed line) and Ej^'^' (dot-dashed line) 
show peculiar behaviors, that consist in alternating slowly decreasing and rapidly decreasing slopes. The average logarithmic 
negativity {E^'^'') (full line) summarizes the behavior of the global entanglement. In panel II we plot the negativities for the 
tauonic neutrino; The behaviors of the negativies for the tauonic instance are similar to the negativities for the muonic instance. 
The curves associated to a given bipartition are plotted with the same plotstyle. The mixing angles Ofj'^^^ and the squared 
mass differences Amfj are fixed at the experimental values (I39p and (|46|l . respectively. We assume the values E = 10 Gel/, 
Op = IGeV, and ^ = for the remaining parameters in Eq. (I45p . The CP-violating phase 5 is put to zero. The x axis is in 
logarithmic scale, and the dimensions are meters. 
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FIG. 7; (Color online) The logarithmic negativities E)^'^' (panel I) and E)^'^' (panel II) as functions of the distance x for 
different choices of the CP- violating phase 5: (a) (5 = (dotted line); (b) 5 = f (dashed line); (b) 5 = n (dot-dashed line). 

iSjj^'^''^' is independent of 5. The mixing angles Ofj^^^ , the squared mass differences Am^^, the parameters E, ap, and ^ are 
fixed as in Figs. [S] and IB] The x axis is in logarithmic scale, and the dimensions are meters. 
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APPENDIX A: ENTROPIC MEASURES FOR THE STATES \W^^\S)) 

Below we report the analytical expressions for the eigenvalues corresponding to the reduced density matrices of 
the states |Vri*^(J)) (a — e,ii,T,s). Let us denote by A^*'-'"'^''^ and A^*'-'''^''-' the eigenvalue vectors associated with the 



16 

reduced density matrices Trj,k,i[\wi^\^)){wi^\^)\] and Trk,i[\wi^\^)){wi^\^)\], respectively. We get 

^(l;2,3,4) ^ ^(2;l,3,4) ^ ^(3;1,2,4) ^ ^4;1,2.3) ^ ^(4;1,2,3) ^ ^(4;1,2,3) ^ ^(4;1,2,3) ^ 1 ^3 ^ ^ (^1) 

^(i;2,3,4) _ 3g |25 — 6cos5i4 — 6cos^23 — 2cos((5i4 + ^23) , H + 6cos5i4 + 6C0S^23 + 2cos(^i4 + (523)}, (A2) 
A(f'i'3'4) = -L|ii-6cos5i4-6cos^23 + 2cos(5i4 + fe),25 + 6cos5i4 + 6cosfe-2cos(^i4 + (523)}, (A3) 
A|f'^'''^) = ^{5-4cos((5i4+(523),31 + 4cos(5i4+(523)}, (A4) 



^U;2,3,4) _ -i- 1 16 - 6 cos 5i4 - 6 COS 623 — 2 cos(5i4 + fe) + 6 cos(5i4 - 534) - 6 cos((5i4 - ^23 - ^34) 

-9 cos 534 + 6 C0S(523 + <^34) + 3 C0S(2^23 + ^3i) ,56 + 6 COS ^14 + 6 COS (523 + 2 C0s(Ji4 + (523) 
-6 COs((5l4 - (534) + 6 COs((5l4 - (523 - (534) + 9 COS ^34 - 6 COs(fe + (534) - 3 C0s(2fe + (534) I , (A5) 



^(2;1,3,4) _ ^jsS — 6c0S(5i4 — 6 COS (523 + 2 COs((5i4 + (523) " 6 C0s(5i4 — i534) — 6 C0s((5i4 — (523 ~ <53 

— 9 COS (534 — 6C0S((523 + (534) + 3 C0S(2(523 + (534) , 16 + 6 c0S(5i4 + 6 COS (523 — 2 c0s((5i4 + (523) 
+6cOs((5i4 - 534) + 6cOs((5i4 - ^23 - <534) + 9 COS (534 + 6C0S((523 + (534) - 3 C0S(2J23 + <534)| , (A6) 

A^.^'^'^'^) = ^{11 + 2cos((5i4 + (523) - 6cos(,5i4 - S^i) - 6cos((523 + ^534) , 

25-2 cos((5i4 + ^23) + 6 cos((5i4 - ^34) + 6 cos(J23 + <534)} , (A7) 

^(i;2,3,4) _ — 156 + 6cos(5i4 + 6 cos (523 + 2cos((5i4 + (523) + 6cos((5i4 — 5m) — 6cos(^i4 — (523 — ^34) 

— 9 COS 634, — 6 C0S(^23 + <534) + 3 C0S(2^23 + 1534) ,16 — 6 COS ^14 — 6 COS (523 ~ 2 C0S(^14 + 823) 

-6 cos((5i4 - ^34) + 6 cos((5i4 - ^23 - S34) + 9 cos ^34 - 6 C0S(J23 + S34) -3 C0S(2^23 + <^34) I , (A8) 

^(2;1,3,4) _ ^|l6 + 6cOS(5i4 + 6C0S(523 - 2cOs((5i4 + 823) ~ 6cOs((5i4 - S34) - 6c0s((5i4 - S23 - S34) 

— 9 COS (534 — 6C0S((523 + S34) + 3 C0S(2(523 + (534) ,56 — 6 COS(5i4 — 6 cos (523 + 2 COs((5i4 + (523) 

+6 cos((5i4 - ^34) + 6 cos((5i4 - ^23 - 634) + 9 cos ^34 + 6 cos( J23 + (534) - 3 cos(2(523 + <534) | , (A9) 



AI^'^'^'^) = ^ {25 - 2 C0S((514 + (523) - 6 C0S(,514 - (534) - 6 C0S((523 + (534) , 

11 + 2 COs((5l4 + fe) + 6 COs((5l4 - ^34) + 6 C0S(J23 + '534)} , (AlO) 



^(1,2;3,4) ^ ^(1,3;2,4) ^ ^(M;2,3) = ^ |0 , Q , 1 , 1} , (All) 

A^''''^'^^ = ^{0,0,7-2cos(5i4 +(523),ll + 2cos(5i4 + (523)}, (A12) 

^(1.3;2,4) ^ -L |o , , 10 - 3 COS (5i4 - 3 cos ^23 + C0s(Ji4 + (523) , 8 + 3 COS Si4 + 3 COS (523 - COs((5i4 + (523) } j(A13) 

^(1.4;2,3) ^ _L |o , , 8 - 3 COS (5l4 - 3 COS (523 - COs(5l4 + 623) ,10 + 3 COS 5l4 + 3 COS + C0S(^14 + (523) } j(A14) 
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^a>2;3,4) ^ _L|o ,0, 10 + COs(Ji4 + S23) - 3cOs((5i4 - S^i) - 3C0S((523 + Ssi) ,8 - COs((5i4 + ^23) 

+ 3 COs((5i4 - ^34) + 3 008(^23 + <534) I , (A15) 

^a:3;2,4) ^ _L|o ,0 , 38 - 6cos(5i4 - 6 COS (523 + 2cos((5i4 + 823) + 6cos((5i4 - ^34) - 6cos((5i4 - S23 - Ssi) 

-9 008^34 - 6C0S((523 + 534) + 3 COs(2(523 + 634) ,34 + 6cOs5i4 + 6cOS(523 - 2cOs((5i4 + 623) 

+6 cos((5i4 - (534) + 6 cos((5i4 - ^23 - (534) + 9 COS (534 + 6 cos(J23 + (534) - 3 cos(2(523 + (534) I , (A16) 

^a>4;2,3) ^ _L|o ,0 , 34 - 6cOs5l4 - 6 COS (523 " 2 C0s(^i4 + (523) + 6 C0s(^i4 - S34) - 6 COs((5i4 - S23 - S34) 

-9 cos (534 + 6cos((523 + (534) + 3cos(2(523 + 834) ,38 + 6cos5i4 + 6 cos (523 + 2cos(Ji4 + 623) 

-6 cos((5i4 ~ (534) + 6 cos((5i4 - 623 - 634) + 9 cos 634 - 6 cos(523 + (534) - 3 cos(2(523 + (534) I , (A17) 

Ai^'^'^'^) = , , 8 ~ C0S((514 + (523) - 3 C0S((514 - (534) - 3 C0S(523 + 534) , 10 + COs((5i4 + (523) 

+ 3C0S((514 - (534) + 3C0S((523 + (534)} , (A18) 

^(1,3;2,4) ^ -L|o ,0 , 34 + 6C0S^14 + 6 COS (523 " 2 C0s(^i4 + S23) - 6 COs((5i4 - (534) - 6cOs((5i4 - (523 ~ ^34) 

—9 COS (534 — 6 cos((523 + ^34) +3 cos(2(523 + (534) , 38 — 6 cos ^14 — 6 cos S23 + 2 cos((5i4 + S23) 

+ 6 COs((5i4 - (534) + 6 COs((5i4 - 623 - 634) + 9 cos 634 + 6 COs(J23 + (534) - 3 C0S(2(523 + (534) I , (A19) 

^(i,4;2,3) ^ -^|o ,0 , 38 + 6cos(5i4 + 6 COS (523 + 2cos((5i4 + S23) + 6cos((5i4 - (534) - 6cos((5i4 - (523 - (^34) 
—9 cos S34 + 6 cos((523 + (534) +3 cos(2(523 + (534) , 34 — 6 cos ^14 — 6 cos S23 — 2 cos((5i4 + S23) 
-6 cos((5i4 - S34) + 6 cos((5i4 - §23 - 634) + 9 cos (534 - 6 cos(J23 + S34) -3 cos(2(523 + (534) I . (A20) 

The von Neumann entropies can be easily written as 

^^^L = -E^"(")log2AiHn). (A21) 
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